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Abstract 

Information about the physical parameters of neutron stars is 
obtained from pulsar observations. The energy balance of the Crab 
Nebula and the Vela X remnant allows one to derive limits for the 
masses of the Crab and Vela pulsars. Glitch observations provide 
further clues on the masses of these two pulsars. The degree of 
confidence with which one should believe the derived numbers is pointed 
out. The possibility to observe neutron stars in binary systems as 
pulsating x-ray sources is discussed. Finally, the importance of 
observing redshifted gamma ray lines from the surface of neutron stars, 
and thus directly measuring either individual or statistical properties 
of these objects is pointed out. 
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I . Introduction 

Besides the well-known chain of arguments that leads in an elimination 
process (Cameron & Maran, 1969) to the generally accepted conclusion that 
pulsars are rotating neutron stars, there is another very gratifying 
result which stems from the comparison of neutron star models and pulsar 
observations. It turns out that for any reasonable equation of state 
the neutron star models, computed as outlined in the preceding talk, 
have moments of inertia, which are just of the right order of magnitude 
to explain the observed energy input into .the Crab Nebula as the loss 
of rotational energy from a rotating neutron star with the parameters of 
the Crab pulsar PSR 0531 + 21. This is a major triumph of the inter- 
action of theory and experiment in this field. It has already been 
mentioned in the talks of H. A. Bethe and D. Pines, but this information really 
belongs to the subject matter of this talk and has therefore been pointed 
out again. We shall discuss in detail the energy balance of the Crab and 
Vela nebulae later. 

Except for this important fact, the interaction between theory 
and observation is very slim to date, so much so that no observer is 
participating in this conference and that observations have to be 
discussed by theorists, one of which has never even looked through a telescope. 
But this is understandable because the wealth of pulsar observations is 
overwhelming, and a comprehensive understanding still has not been 
achieved. Only the puzzling flitches* of the Crab and Vela pulsar have 
been well incorporated into various competing theories. These will be 
discussed at the end of this talk. 

Preceding page blank 
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In the following we shall attempt to gather some of the pieces 
in the puzzle that links neutron star models and pulsar 

observations. We shall try very ambitiously to derive 

parameters for the rotating neutron stars that are represented by the 
names Crab and Vela in the astronomical observations. 

II. Radio Pulsars 

All of the 61 pulsars listed by Manchester & Taylor (1972), except 
two (Crab and Vela), are seen only by their pulsed radio emission. It 
is quite obvious that the pattern, shape and polarization of these radio 
pulses contains quite a lot of information on the internal structure of 
the pulsar. We are, however, not yet able to understand the message we 
are getting. No satisfying quantitative description of the electro- 
magnetic link between the rotating neutron star and the radiation pattern 
of the pulsar has been given so far. Indeed, not even the case of a 
magnetosphere of radiating particles, where the axis of the magnetic field 
coincides with the rotation axis of the neutron star, has been solved. 
Whereas to explain the pulse producing mechanism one would have to treat 
the much more complicated case of at least a slight deviation from axial 
symmetry . 

In the absence of a detailed mechanism for pulsar emission, only energy 
considerations can be employed to obtain information on the physical 

properties of the rotating neutron star. For 22 of the 61 listed pulsars 
(Manchester 6c Taylor, 1972) both frequency 0 and change of frequency 
H have been measured. Then by determining their rate of loss 

of energy 


i no 
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we could in principle find the moment of inertia I of these neutron 
stars. This in turn would precisely fix the mass and density profile of 
the star once an equation of state has been chosen. 

The energy in radio pulses E pu ^ g probably is only a small part 
of the total energy release, and therefore only very crude limits on the 
physical parameters of a neutron star may be derived. Even this modest 
undertaking does not look very promising, however, since neither has the 
flux in radio pulses , a well-defined mean value for a specific pulsar, 
nor is the distance to pulsars known to a high accuracy. 

The flux in radio pulses is known to fluctuate very strongly on a 
scale of weeks. This makes it very difficult to define a quantity like 
the mean pulse intensity at a given frequency for any pulsar. Furthermore 
the distance to pulsars can only be estimated (again excepting Crab and 
Vela) by their dispersion measure, which gives the average value of the 
electron density along the line of sight. To derive the distances of 
pulsars from the dispersion measure one would need a precise knowledge 
of the interstellar medium. In reality, however, we use the pulsars to 
obtain more information on the interstellar medium; and, it has been 
found that in a region of 100 pc around the sun, the average electron 

— o 

density is almost 0.1 cm . This is roughly comparable 

to the average density of the atomic hydrogen in that region, and 
indicates that we are surrounded neither by a classical HII region nor 
by an HI region (c.f. Biermann, 1972). 

Thus there is virtually no feedback to the neutron star models 
through the observations of radio pulses. If we nevertheless 
make some crude estimates of the energy in radio pulses for various 
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pulsars, we find Ep U ^ s ^ 10 — ergs/sec. While that result is laden 

will all the uncertainties discussed above, it permits one rather 
vague, but still interesting conclusion: assuming moments of inertia 

between 10^ and 10^ g cm^ we can compute the rotational energy output 

E ro t = I ^ for various pulsars. We then find that the energy in 

• -9 * 

radio pulses for the Crab pulsar Ep U ]^ s ^ 10 ^ ro t > w ^ll e f° r many of 
the older pulsars Ep U -^ s is a much bigger fraction of E rofc , such as 
PSR 0809 + 74, where Ep U ^ g ^ E rot . This seems to indicate that as 
pulsars grow older they spend a bigger and bigger fraction of their 
rotational energy output in the production of radio pulses. 

Ill . The Crab Pulsar 

The Crab pulsar PSR 0531 + 21 is located at the center of the Crab 
Nebula, the remnant of a supernova that exploded in 1054. Because of 
its location in the nebula a distance estimate independent of the 
dispersion measure can be obtained, and furthermore the energy balance 
of the nebula can be used to get more information on the energy output 
of this pulsar. In addition, pulses from this object have been observed 
not only in the radio regime but also in the optical and x-ray frequencies. 

The supernova of 1054 is widely considered to have been of type I, 
but Minkowski (1968) has case doubt on that. Trimble (1968) has compared 
the radial velocities and proper motions in the thick filamentary shell, 
which is a projected elliptical object with semiaxes of 3* and 2 1 , and 
obtained a distance of 2 kpc on the assumption that the 3-dimensional 
nebula is a prolate ellipsoid. Trimple & Woltjer (1971) have pointed 
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out that the uncertainties in this value are large and that a distance 
as low as 1.2 kpc or as high as 2.5 kpc cannot be excluded. 

Let us, however, be definite and adopt a distance of 2 kpc 

for the Crab Nebula. This value seems to be the most widely used. 
Observations of the nebula indicate that energy must be supplied to it 
continuously. Assuming that the Crab pulsar is the only source of 
energy in the nebula, one can determine limits on the energy output 
of the pulsar (Rees & Trimble, 1970; Borner & Cohen, 1972) by considering 
the energy balance. 

The only well established energy loss is the synchrotron radiation 

• Oo 

which implies E S y nc ^ “ 1.2 x 10 ° erg/sec (Baldwin, 1971), if the distance 
to the Crab is 2 kpc. The pulsar has to replenish at least 

the electrons producing the optical and x-ray synchrotron radiation, 
because these particles have half-lives of less than 100 years. So a 
rough estimate obtained from the observed spectrum (Baldwin, 1971) 

OQ 

indicates that the pulsar has to supply continuously at least 0.8 x 10 J 
erg/sec. For the pulsar to replenish this energy, via loss of 
rotational energy 


E rot = inn > 

its moment of inertia has to be at least 1.8 x 10^ g cm2 (c.f. Cohen & 
Cameron, 1971), corresponding to line (a) in Figure 7 of the preceding 
talk. Neutron star models corresponding to various equations of state 
have been discussed in the preceding talk and it can be seen from 
Figure 2 that a model with 1 = 1.8 x 10^ g cm^ has a mass of 

0.34 Ms (BPS), 0*. 36 Mo (BJ), 0.36 Mo (BBS), 0.26 Mo (CCLR) . 
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There BPS corresponds to the equation of state published in Baym et al. 

(1971); BJ: Bethe & Johnson (1973); BBS: Bethe et al. (1970); CCLR: 

Cohen et al. (1970). All the equations of state discussed in the 
preceding talk, except one, can easily provide neutron star models big 
enough to exceed this lower limit. The exception is the work of Leung 6c Wang 
(1971); the maximum moment of inertia in their equation of state, numbered (I) 

A/l 2 

is less than 0.2 x 10r^ g cm . That is the Crab pulsar is definitely not 
among the stable neutron stars they compute using their equation of state, 
which does not incorporate repulsion between baryons. Even equation (II) 
of Leung & Wang (1971), which gives I max = 1.05 x 10^ g cm2 i s too low 
(Leung 6c Wang in Figure 7 of preceding talk), although some repulsion is 
assumed to be present in this case (II). This indicates, as has been 
discussed already from the nuclear physics point of view by H. A. Bethe, 
that the repulsion between the nucleons and hyperons at short distances 
plays an important role which may be the dominant feature at high densities. 

It is interesting to see that this is also suggested by astrophysical evidence. 

The evidence in favor of repulsive interactions becomes even stronger, 
if we take into account the protons that are pulled from the surface of the 
rotating neutron star by strong electric fields in the model proposed by 
Goldreich 6c Julian (1969). According to this model these protons are 
accelerated to the same energies as the electrons producing synchrotron 
radiation. In this case the minimum energy loss is twice the synchrotron 
radiation loss. This has the consequence that the minimum moment of inertia 
of the Crab increases by a factor of 2: 

^in > 4 x 1q44 § cm2 

The neutron star mass is then > 0.5 (BBS, BJ). Thus if this model is 

valid we obtain the result that the Crab pulsar is a neutron star with a 
mass of at least 0.5 Mq. 
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The acceleration mechanism of Gunn & Ostriker (1969) can also be 
tested. They require that the protons get ten times the energy of the 
electrons. Assuming the fluxes of electrons and ions from the pulsar 
to be equal, this would require ten times the synchrotron energy for 
the protons, leading to a moment of inertia of 2 x 10 ^ g C m^ (line (c) 
in Figure 7 of the preceding talk). Only the CCLR equation of state has 
models with moments of inertia of that magnitude. But even then this 
condition is satisfied only over a small density range near the mass peak. 
It therefore seems that the model of Gunn & Ostriker (1969) should be 
modified quanitatively . 

Conclusions based on energy losses from the Crab nebula due to 
the expanding supernova shell (Rees & Trimble, 1970; Borner 6c Cohen, 

1971) are much more uncertain than the preceding considerations. 
Observations of the filaments in the expanding supernova shell 
(Woltjer, 1958; Trimble, 1968) show that the expansion velocity at 
present is higher than would correspond to an expansion at constant 
velocity since 1054. It seems that the nebula is accelerating now with 
an acceleration of 

v = 0.0014 cm/sec^ . 

The nebula might, however, be decelerating now with the velocity still 
higher than the average, if it had been accelerating rather strongly in 
the past. The energy of the expanding supernova shell would change due 
to this acceleration at a rate 

E acc = J P v v d(Vol) (1) 

p: density in the nebula. 

It would also change by the 'snow plow' effect, the change in mass of the 
supernova shell as interstellar material piles up along the rim 
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E 1 =1/2 Pm v A 

plow 


A: surface area of nebula 


p^: density of interstellar material. 


Since we do not know whether the supernova shell is accelerating or 

decelerating at present, we investigate both cases. If the shell is 

decelerating, the energy gained by deceleration will be spent in the 

f snow plow 1 effect described by (2), and perhaps totally balance it. 

Thus E , . n =E +E.. = 0 is a distinct possibility. No 

shell acc plow 

further limits on the parameters of the Crab pulsar except those derived 
earlier from synchrotron radiation can be found in this case. We should 
notice, however, that in principle one could directly measure the value 
of v at present, and thus decide the question of acceleration or 
deceleration. If the currently accepted values for acceleration and 
snow-plow are used, we find 

• OQ 

E acc = 1.6 x 10 J erg/sec (3) 


"plow = 1.7 X 10 8 erg/sec 


hence E 


shell 


= 3.3 x lCr° erg/s'ec 


This energy has to be supplied either directly by the pulsar via the 
low frequency waves emitted or by the adiabatic expansion of a 
relativistic gas (Trimble & Rees, 1970). In both cases a rotating 
neutron star with a moment of inertia of 10 g cm^ can continuously 
supply that energy (line (b) in Figure 7 of the preceding talk). The 


equations of state that can furnish a neutron star model with a moment of inertia 


big enough (BBS, BJ, CCER) to give a mass of 1.2 Mq (BBS, BJ) for this model. 
If little material was lost during the collapse, the star would have 
had a mass of 1.35 prior to the collapse, above the Chandrasekhar 
limit for typical white dwarfs. It should be remembered that these 



10 


energy losses are rather uncertain, and that the pulsar has to supply the 
energy continuously only if the energy content of the gas of relativistic 
particles in the nebula is maintained at its present level. All the 
uncertainties can, however, be decided by future observations, and 
thus a value of 1.2 Mq for the Crab pulsar may be confirmed some day 
with a much higher degree of confidence than we have now. 

If the shell is pushed out by low-frequency waves from the pulsar, 
and if according to Gunn & Ostriker (1969) protons get ten times the 
energy of the electrons, then only the maximum mass (near the mass peak) 
neutron star models of the CCLR equation of state can fulfill this requirement, 
as indicated by line (d) in Figure 7 of the preceding talk. This particle 
acceleration mechanism therefore seems lo be unrealistic. 

If the distance to the Crab were less or more than 2 kpc, the 
limits derived above would have to be scaled down or up accordingly. 

If we go to the extreme values for the distance of 1.2 kpc to 2.5 kpc, 
the moment of inertia necessary to account for the short-lived 

o 

synchrotron particles would vary between 0.8 and 2.4 x 10 g cm^ 

(0.2 Mq and 0.4 M 0 , respectively (BBS)). 

A final remark might be of interest. It pertains to the suggestion 
that cosmic rays are produced in the electromagnetic field of pulsars. 

The spectrum of galactic cosmic rays above 300 MeV shows about 50 times 
as much energy is present in protons than in electrons. It is clear 
from the foregoing discussion that the Crab pulsar could not produce 
such a ratio of proton to electron array. This throws considerable doubt 
on the hypothesis that the high energy galactic cosmic rays are all 


produced by pulsars. 
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IV. The Vela Pulsar 

The Vela pulsar PSR 0833-45 is associated with the supernova 

4 

remnant Vela X, which is about 1.1 x 10 years old, has a radius of 
10 pc, and is at a distance of 500 pc (Milne, 1970). 

If the electromagnetic radiation emitted by Vela X is synchrotron, 
then from synchrotron theory the energy content in the gas can be 
estimated at ~ 10 ^ (Tucker, 1971). Assuming constant velocity, if the 
expansion velocity of the supernova shell were constant, which would mean 
v = 880 km/sec, the energy loss through adiabatic expansion would be 
E a d = 2.4 x 1037 erg/sec (6) 

On the other hand the rotational energy lost by a neutron star with 
the parameters of PSR 0833-45 is between 

E rot = 4 x 10^6 erg/sec to 8.5 x 10^ erg/sec (7) 

for models with mass between 0.8 Mq and 1.7 Mq (BBS). Although this 
is rather large compared to the loss of 4 x 10^5 erg/sec in x-rays 
(Tucker, 1971) and 10^3 erg/sec in radio (Milne, 1970), the pulsar 
cannot supply the energy given in equation (6), The shell must 
therefore have been decelerating. Let us assume that the deceleration 
at present is very small, and that the main contribution to E s hell is 
by 'snow plow 1 (Borner 6c Cohen, 1972). If it is further assumed that a 
neutron star of 1.2 Mq is present to balance the expansion losses, then 
the velocity of expansion can be determined to be 240 km/sec (Borner 6c 
Cohen, 1972). It is amusing to note that at about the same time Wallerstein 
6c Silk (1971) independently (neither group knew of the other's work until 
after publication) measured the expansion velocity of Vela X by an observation 
of Call lines in that direction, and they found precisely that value of 
240 km/sec for the expansion velocity. 
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V. Pulsar Glitches 

The speed-ups of both Crab and Vela have been discussed by D. Pines, 
along with the various theoretical attempts to account for them. The 
’starquake’ theory advocated by Pines in his talk, as well as the 
’accretion 1 model (Borner & Cohen, 1971) are the only theories of 
pulsar glitches that derive limits on the mass (and other parameters) 
of Crab and Vela. Both theories use the two-component model for the 
interaction between the crust-charged particle system and the neutron 
superfluid to account for the relaxation phenomena of the post-speed-up 
behavior of the pulsar. They differ in the mechanism evoked to produce 
the initial glitch. Both theories are not without problems. 

The starquake theory as described in Pines’ talk views the smaller 
glitches of the Crab pulsar ~ 10“^ j as a sudden relaxation of elastic 
stresses which accumulate in the rigid outer layers, the ’crust', formed 
by a lattice of nuclei, of the pulsar as it is slowing down. The 
magnitude of this effect depends on the magnitude of the stresses that 
can be built up in the crust. Baym & Pines (1971) used an equation of 
state developed by Baym, Bethe & Pethick (1972), where nuclei in the 
lattice become very large, up to Z ~ 200. On the other hand, Negele, 
in his talk, came to the conclusion that Z ~ 40 is to be expected in 
the crust of neutron stars. The amount of stress that can be built up 
increases monotonically with 2. It therefore seems that 
the crustquake theory has serious problems in explaining even 
the Crab pulsar glitches, if the neutron star crust consists only of 
small nuclei. Even for Z ~ 200 nuclei in the lattice, the requirement 
to have a typical Crab speed-up every two years leads to a picture of 
the Crab pulsar as an almost completely solidified star. The crustquake 
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theory predicts a mass of the Crab of less than 0.15 This limit 

should not be taken too much at face value, but it illustrates the 
difficulties of this theory. A neutron star of 0.15 Hq might not 
even be formed in a supernova explosion (c.f. preceding talk), because 
its binding energy is so low that it becomes energetically more 
favorable to form dispersed Fe56. Furthermore 0.15 disagrees with 
all the observational limits discussed above. The larger glitches of 
the Vela pulsar ^ ~ 10"^ j are ascribed to corequakes by Pines, sudden 

relaxations of stress stored in the neutron star's central core made of 
a hadron lattice. The question of whether or not a solid hadron core 
may exist in neutron stars has been discussed extensively at this meeting, 
but no definite conclusion has been reached. The work of Canuto & 


Chitre (Canuto ! s talk at this meeting) suggests that a solid core 

might form at densities above 1.5 x 10^ g/cm^. The corequake 

theory therefore predicts that the Vela pulsar has a central density 

of at least 1.5 x 10^ g/cm^, which makes Vela a rather heavy neutron 

star with mass greater than 0.8 (BPS), 1.2 (BBS; CCLR), and 1.5 Mq (BJ) 

according to the various equations of state. The stress in the core is 

not built up between glitches in this model, but rather each glitch 

takes out a small part of a huge reservoir of elastic stress stored in 

the core. One has to explain why this produces glitches of ( ^2 ~ 10“^), 

\ n / 

instead of a continuous relaxation or one extremely big jump. 

In the accretion model (Borner 6c Cohen, 1971) the initial glitch 

'aq\ 

— ] is ascribed to the infall of material onto the neutron star, 

W J G 

transferring angular momentum to the crust and speeding it up. After 
some time the initial speed-up of the crust is transferred to the interior, 
and the pulsar settles down to a long-term frequency increase : 
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= Aj r 


Here we have J c : angular momentum of crust 


moment of Inertia of crust 


AQ Aj f 


By choosing a definite neutron star model, everything is determined 
from the measured quantities ( and (— \ . Even the infalling mass 


from the measured quantities [ — and — \ . Even the infalling mass 

\ n /g \n 

10 

Am can be found. In this simple model we find Am ~ 10" for crab 

glitches, Am ~ 10" Mq for Vela speed-ups. It is quite clear that 
initial conditions can be formulated which would exactly produce the 
observed behavior. A massive body flung off from the vicinity of the 
pulsar at some early stage in its life, not quite reaching escape 
velocity, but making just one loop and falling back onto a slowed down 
neutron star, would certainly transmit the right amount of angular 
momentum. It would also be able to accrete easily on the neutron star 
surface because it would have little excess angular momentum. The mass 
balance certainly is no problem, because Am is very small compared to 
typical pulsar masses. The big question is, however, whether the 
conditions prevailing after a supernova explosion can lead (with a certain 
non-zero probability) to the initial conditions needed for the accretion 
model. There certainly is a lot of homework to do, but not withstanding 
these theoretical difficulties we may point out that the accretion model 
agrees well with the observations. And it uses the same mechanism to 
explain the glitches of both Crab and Vela pulsars. 

Viewing and as quantities determined (with a considerable 

uncertainty) by the observations, one derives from (8) and (9) the 


condition 
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( 10 ) 


This leads to the condition that the mass of the Crab pulsar be greater 

than 1 Mq; for Vela (10) does not impose any restrictions on the neutron 
star mass . 

VI. Other Possibilities to Observe Neutron Stars 



Recent observations from Uhuru have established the binary nature 

of two periodic pulsating x-ray sources: Cen X-3 (Schreier et al., 1972) 

with a period of 4.84 sec and Her X-l (Tananbaum et al., 1972) with a 

period of 1.24 sec. A model may be suggested in which the x-ray source 

is a neutron star emitting x-rays by radiating as a black body in a 

number of hot spots. The rotation period of the neutron star provides 

the timing mechanism of the pulsation, and the x-rays are produced by 

the accretion of mass from a binary companion. Neutron stars accreting 

matter as models for x-ray sources have been proposed already by a 

number of authors (e.g., Shklovsky, 1967), but a reexamination of the 

older proposals in the light of these new observations might be worthwhile. 

We will not deal with the many intricate questions involved in such a 

model (c.f . Borner et al., 1972), but just assume it to be valid and use 

M2 3 o 

the astronomical observation of the mass function sin- 3 i 

(Mi + M2) 2 

to get information on the mass of the neutron star. An optical identification 
of the main star will give its mass within rather narrow limits. The 
inclination sin i may be determined from an analysis of the pulse shape 
(Borner et al. (1972) find sin i w 1 for Her X-l). Then the mass 
function will directly determine the mass M]_ of the neutron star. If 
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e.g. for Her X-l, where ^ + M 2) 2 ~ 0,85 ^ for Sin 1 ~ M 2 is an 
F-type star of ~ 2 Mq, then it follows ^ 1 Mq . 

(VI. 2) Gamma Ray Lines from Old Neutron Stars 

Johnson et al. (1972) have reported low energy gamma ray observations 
from the galactic center region, showing a statistically significant 
spectral feature at 473 + 30 keV. This line emission with a total 
photon flux of 1.8 x 10*3 cm“2 sec’l has been interpreted by Ramaty 
et al. (1972) as gravitationally redshifted positron annihilation 
radiation from the surface of old neutron stars. The production of 
positrons is attributed to nuclear reactions on the neutron star surface, 
induced by the accretion of interstellar material. It is found that for an 
accretion rate of 10-*- 1 g/sec, or 6 x 10 3 ^ particles/sec, the redshifted 
positron annihilation yield of a single neutron star is 1.2 x 1C)33 photons/ 
sec. Thus, in order to account for the observations, a total of 1.5 x 10^® 
old neutron stars in the galactic center region is required. Positron 
annihilation radiation is normally at 511 keV, so this interpretation 
requires redshifts ranging from 0.016 to 0.13. Since a given redshift 
completely determines a specific neutron star model, we have as a 
consequence that the majority of these old neutron stars have masses of 
less than 0.8 Mq . While the numbers are probably quite uncertain still, 
there appears the interesting possibility that observations of redshifted gamma 
ray lines provide a direct measurement of the general distribution of 
the physical parameters of neutron stars. The principal observational 
tests of the model of Ramaty et al. (1972) would be the detection of 
nuclear gamma ray lines from the galactic center and redshifted positron 
annihilation radiation from the galactic disk. 
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An accretion of 10 g/sec (assuming a surface composition of CNO 
or heavier nuclei) should produce about l(p9 pos itrons /sec from a single 
neutron star, and (if the object is at a distance of 3 kpc) a flux of 
redshifted positron annihilation radiation at earth of 10“^ protons/cm2 sec 
This gamma ray flux is below the presently available detector sensitivities 
but in the future we might perhaps be able to see the redshifted positron 
annihilation radiation from binary x-ray sources or from single nearby 


neutron stars. 
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APPENDIX 


We reproduce here the Figures 2 and 7 of the preceding talk 
’’Hadron Star Models” by Jeffrey M. Cohen and Gerhard Borner. 

Fig 2. Gravitational mass vs. density for various equations of state. 
Fig. 7. Moment of inertia (curved lines) of neutron star models as a 
function of density for various equations of state. 
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MOMENT OF INERTIA I VS. CENTRAL DENSITY 






